Multiscale Imaging of Evapotranspiration

Evapotranspiration (ET; evaporation + transpiration) is central to a wide range of biological, chemical, and physical processes in the Earth system. Accurate remote sensing of ET is challenging due to the interrelated and generally scale dependent nature of the physical factors which contribute to the process. The evaporation of water from porous media like sands and soils is an important subset of the complete ET problem. Chapter 1 presents a laboratory investigation into this question, examining the effects of grain size and composition on the evolution of drying sands. The effects of composition are found to be 2-5x greater than the effects of grain size, indicating that differences in heating caused by differences in reflectance may dominate hydrologic differences caused by grain size variation. In order to relate the results of Chapter 1 to the satellite image archive, however, the question of information loss between hyperspectral (measurements at 100s of wavelength intervals) laboratory measurements and multispectral (≤ 12 wavelength intervals) satellite images must be addressed. Chapter 2 focuses on this question as applied to substrate materials such as sediment, soil, rock, and non-photosynthetic vegetation. The results indicate that the continuum that is resolved by multispectral sensors is sufficient to resolve the gradient between sand-rich and clay-rich soils, and that this gradient is also a dominant feature in hyperspectral mixing spaces where the actual absorptions can be resolved. Multispectral measurements can be converted to biogeophysically relevant quantities using spectral mixture analysis (SMA). However, retrospective multitemporal analysis first requires cross-sensor calibration of the mixture model. Chapter 3 presents this calibration, allowing multispectral image data to be used interchangeably throughout the Landsat 4-8 archive. In addition, a theoretical explanation is advanced for the observed superior scaling properties of SMA-derived fraction images over spectral indices. The physical quantities estimated by the spectral mixture model are then compared to simultaneously imaged surface temperature, as well as to the derived parameters of ET Fraction and Moisture Availability. SMA-derived vegetation abundance is found to produce substantially more informative ET maps, and SMA-derived substrate fraction is found to yield a surprisingly strong linear relationship with surface temperature. These results provide context for agricultural applications. Chapter 5 investigates the question of mapping and monitoring rice agricultural using optical and thermal satellite image time series. Thermal image time series are found to produce more accurate maps of rice presence/absence, but optical image time series are found to produce more accurate maps of rice crop timing. Chapter 6 takes a more global approach, investigating the spatial structure of agricultural networks for a diverse set of landscapes. Surprisingly consistent scaling relations are found. These relations are assessed in the context of a network-based approach to land cover analysis, with potential implications for the scale dependence of ET estimates. In sum, this thesis present a novel approach to improving ET estimation based on a synthesis of complementary laboratory measurements, satellite image analysis, and field observations. Alone, each of these independent sources of information provides novel insights. Viewed together, these insights form the basis of a more accurate and complete geophysical understanding of the ET phenomenon

Humans on Earth: Global extents of anthropogenic land cover from remote sensing

This review provides a perspective of the current state of the art in remote sensing of anthropogenic land cover and human-modified landscapes at global scales. The fact that humans have adapted to almost all of Earth’s environments, yet remain strongly clustered within each of these environments influences both the nature of anthropogenic impact on Earth’s landscapes and the challenges of mapping it. Remote sensing provides a consistent synoptic view of these environments by mapping the land cover associated with the anthropogenic land uses of settlement and food production, as well as their complement in forest cover. We give brief descriptions and illustrative comparisons of several current land cover products representing the global extents of settlements, agriculture and forests derived from remote sensing. To accommodate the challenges inherent to mapping any land cover at widely varying scales, we compare size distributions of spatially contiguous land cover (rather than total area) for several global land cover products. Despite the use of different sensors and different mapping criteria, there is remarkable consistency in the size distributions of these products – both within and across land cover class. Rank-size distributions of settlements, agricultural areas and forests are all well-described by power laws spanning more than four orders of magnitude in both area and number. This consistency in the form of the distributions suggests fundamental similarities among different types of land cover. The observed similarities can be explained by depicting land cover mosaics as co-evolving spatial networks sharing common processes of nucleation, growth and connection

Which Vegetation Index? Benchmarking Multispectral Metrics to Hyperspectral Mixture Models in Diverse Cropland

The monitoring of agronomic parameters like biomass, water stress, and plant health can benefit from synergistic use of all available remotely sensed information. Multispectral imagery has been used for this purpose for decades, largely with vegetation indices (VIs). Many multispectral VIs exist, typically relying on a single feature—the spectral red edge—for information. Where hyperspectral imagery is available, spectral mixture models can use the full VSWIR spectrum to yield further insight, simultaneously estimating area fractions of multiple materials within mixed pixels. Here we investigate the relationships between VIs and mixture models by comparing hyperspectral endmember fractions to six common multispectral VIs in California’s diverse crops and soils. In so doing, we isolate spectral effects from sensor- and acquisition-specific variability associated with atmosphere, illumination, and view geometry. Specifically, we compare: (1) fractional area of photosynthetic vegetation (Fv) from 64,000,000 3–5 m resolution AVIRIS-ng reflectance spectra; and (2) six popular VIs (NDVI, NIRv, EVI, EVI2, SR, DVI) computed from simulated Planet SuperDove reflectance spectra derived from the AVIRIS-ng spectra. Hyperspectral Fv and multispectral VIs are compared using both parametric (Pearson correlation, ρ) and nonparametric (Mutual Information, MI) metrics. Four VIs (NIRv, DVI, EVI, EVI2) showed strong linear relationships with Fv (ρ > 0.94; MI > 1.2). NIRv and DVI showed strong interrelation (ρ > 0.99, MI > 2.4), but deviated from a 1:1 correspondence with Fv. EVI and EVI2 were strongly interrelated (ρ > 0.99, MI > 2.3) and more closely approximated a 1:1 relationship with Fv. In contrast, NDVI and SR showed a weaker, nonlinear, heteroskedastic relation to Fv (ρ < 0.84, MI = 0.69). NDVI exhibited both especially severe sensitivity to unvegetated background (–0.05 < NDVI < +0.6) and saturation (0.2 < Fv < 0.8 for NDVI = 0.7). The self-consistent atmospheric correction, radiometry, and sun-sensor geometry allows this simulation approach to be further applied to indices, sensors, and landscapes worldwide

Land Cover Dynamics on the Lower Ganges–Brahmaputra Delta: Agriculture–Aquaculture Transitions, 1972–2017

Aquaculture in tropical and subtropical developing countries has expanded in recent years. This practice is controversial due to its potential for serious economic, food security, and environmental impacts—especially for intensive operations in and near mangrove ecosystems, where many shrimp species spawn. While considerable effort has been directed toward understanding aquaculture impacts, maps of spatial extent and multi-decade spatiotemporal dynamics remain sparse. This is in part because aquaculture ponds (ghers) can be challenging to distinguish from other shallow water targets on the basis of water-leaving radiance alone. Here, we focus on the Lower Ganges–Brahmaputra Delta (GBD), one of the most expansive areas of recent aquaculture growth on Earth and adjacent to the Sundarbans mangrove forest, a biodiversity hotspot. We use a combination of MODIS 16-day EVI composites and 45 years (1972–2017) of Landsat observations to characterize dominant spatiotemporal patterns in the vegetation phenology of the area, identify consistent seasonal optical differences between flooded ghers and other land uses, and quantify the multi-decade expansion of standing water bodies. Considerable non-uniqueness exists in the spectral signature of ghers on the GBD, propagating into uncertainty in estimates of spatial extent. We implement three progressive decision boundaries to explicitly quantify this uncertainty and provide liberal, moderate, and conservative estimates of flooded gher extent on three different spatial scales. Using multiple extents and multiple thresholds, we quantify the size distribution of contiguous regions of flooded gher extent at ten-year intervals. The moderate threshold shows standing water area within Bangladeshi polders to have expanded from less than 300 km2 in 1990 to over 1400 km2 in 2015. At all three scales investigated, the size distribution of standing water bodies is increasingly dominated by larger, more interconnected networks of flooded areas associated with aquaculture. Much of this expansion has occurred in immediate proximity to the Bangladeshi Sundarbans

Spatiotemporal evolution of COVID‐19 infection and detection within night light networks: comparative analysis of USA and China

The spatial distribution of population affects disease transmission, especially when shelter in place orders restrict mobility for a large fraction of the population. The spatial network structure of settlements therefore imposes a fundamental constraint on the spatial distribution of the population through which a communicable disease can spread. In this analysis we use the spatial network structure of lighted development as a proxy for the distribution of ambient population to compare the spatiotemporal evolution of COVID-19 confirmed cases in the USA and China. The Visible Infrared Imaging Radiometer Suite (VIIRS) Day/Night Band sensor on the NASA/NOAA Suomi satellite has been imaging night light at ~700 m resolution globally since 2012. Comparisons with sub-kilometer resolution census observations in different countries across different levels of development indicate that night light luminance scales with population density over ~ 3 orders of magnitude. However, VIIRS’ constant ~ 700 m resolution can provide a more detailed representation of population distribution in peri-urban and rural areas where aggregated census blocks lack comparable spatial detail. By varying the low luminance threshold of VIIRS-derived night light, we depict spatial networks of lighted development of varying degrees of connectivity within which populations are distributed. The resulting size distributions of spatial network components (connected clusters of nodes) vary with degree of connectivity, but maintain consistent scaling over a wide range (5 × to 10 × in area & number) of network sizes. At continental scales, spatial network rank-size distributions obtained from VIIRS night light brightness are well-described by power laws with exponents near −2 (slopes near −1) for a wide range of low luminance thresholds. The largest components (10e4 to 10e5 sq.km) represent spatially contiguous agglomerations of urban, suburban and periurban development, while the smallest components represent isolated rural settlements. Projecting county and city-level numbers of con rmed cases of COVID-19 for the USA and China (respectively) onto the corresponding spatial networks of lighted development allows the spatiotemporal evolution of the epidemic (infection and detection) to be quantified as propagation within networks of varying connectivity. Results for China show rapid nucleation and diffusion in January 2020 followed by rapid decreases in new cases in February. While most of the largest cities in China showed new confirmed cases approaching zero before the end of February, most of these cities also showed distinct second waves of cases in March or April. Whereas new cases in Wuhan did not approach zero until mid-March, as of December 2020 it has not yet experienced a second wave of cases. In contrast, the results for the USA show a wide range of trajectories, with an abrupt transition from slow increases in confirmed cases in a small number of network components in January and February, to rapid geographic dispersion to a larger number of components shortly before mobility reductions occurred in March. Results indicate that while most of the upper tail of the network had been exposed by the end of March, the lower tail of the component size distribution has only shown steep increases since mid-June

Global cross-calibration of Landsat spectral mixture models

Data continuity for the Landsat program relies on accurate cross-calibration among sensors. The Landsat 8 Operational Land Imager (OLI) has been shown to exhibit superior performance to the sensors on Landsats 4–7 with respect to radiometric calibration, signal to noise, and geolocation. However, improvements to the positioning of the spectral response functions on the OLI have resulted in known biases for commonly used spectral indices because the new band responses integrate absorption features differently from previous Landsat sensors. The objective of this analysis is to quantify the impact of these changes on linear spectral mixture models that use imagery collected by different Landsat sensors. The 2013 underflight of Landsat 7 and Landsat 8 provides an opportunity to cross calibrate the spectral mixing spaces of the ETM+ and OLI sensors using near-simultaneous acquisitions of radiance measurements from a wide variety of land cover types worldwide. We use 80,910,343 pairs of OLI and ETM+ spectra to characterize the Landsat 8 OLI spectral mixing space and perform a cross-calibration with Landsat 7 ETM+. This new global collection of Landsat spectra spans a greater spectral diversity than those used in prior studies and the resulting Substrate, Vegetation, and Dark (SVD) spectral endmembers (EMs) supplant prior global Landsat EMs. We find only minor (− 0.01 < μ < 0.01) differences between SVD fractions for coregistered pairs of spectra unmixed using the new sensor-specific endmembers identified in this analysis. Root mean square (RMS) misfit fractions are also shown to be small (< 98% of pixels with < 5% RMS), in accord with previous studies using standardized global endmembers. Finally, vegetation is used as an example to illustrate the empirical and theoretical relationship between commonly used spectral indices and subpixel fractions. We include the new global ETM+ and OLI EMs as Supplementary Materials. SVD fractions unmixed using global EMs thus provide easily computable, linearly scalable, physically based measures of subpixel land cover area which can be compared accurately across the entire Landsat 4–8 archive without introducing any additional cross-sensor corrections

Spectral Mixture Analysis as a Unified Framework for the Remote Sensing of Evapotranspiration

This study illustrates a unified, physically-based framework for mapping landscape parameters of evapotranspiration (ET) using spectral mixture analysis (SMA). The framework integrates two widely used approaches by relating radiometric surface temperature to subpixel fractions of substrate (S), vegetation (V), and dark (D) spectral endmembers (EMs). Spatial and temporal variations in these spectral endmember fractions reflect process-driven variations in soil moisture, vegetation phenology, and illumination. Using all available Landsat 8 scenes from the peak growing season in the agriculturally diverse Sacramento Valley of northern California, we characterize the spatiotemporal relationships between each of the S, V, D land cover fractions and apparent brightness temperature (T) using bivariate distributions in the ET parameter spaces. The dark fraction scales inversely with shortwave broadband albedo (ρ < −0.98), and show a multilinear relationship to T. Substrate fraction estimates show a consistent (ρ ≈ 0.7 to 0.9) linear relationship to T. The vegetation fraction showed the expected triangular relationship to T. However, the bivariate distribution of V and T shows more distinct clustering than the distributions of Normalized Difference Vegetation Index (NDVI)-based proxies and T. Following the Triangle Method, the V fraction is used with T to compute the spatial maps of the ET fraction (EF; the ratio of the actual total ET to the net radiation) and moisture availability (Mo; the ratio of the actual soil surface evaporation to potential ET at the soil surface). EF and Mo estimates derived from the V fraction distinguish among rice growth stages, and between rice and non-rice agriculture, more clearly than those derived from transformed NDVI proxies. Met station-based reference ET & soil temperatures also track vegetation fraction-based estimates of EF & Mo more closely than do NDVI-based estimates of EF & Mo. The proposed approach using S, V, D land cover fractions in conjunction with T (SVD+T) provides a physically-based conceptual framework that unifies two widely-used approaches by simultaneously mapping the effects of albedo and vegetation abundance on the surface temperature field. The additional information provided by the third (Substrate) fraction suggests a potential avenue for ET model improvement by providing an explicit observational constraint on the exposed soil fraction and its moisture-modulated brightness. The structures of the T, EF & Mo vs SVD feature spaces are complementary and that can be interpreted in the context of physical variables that scale linearly and that can be represented directly in process models. Using the structure of the feature spaces to represent the spatiotemporal trajectory of crop phenology is possible in agricultural settings, because variations in the timing of planting and irrigation result in continuous trajectories in the physical parameter spaces that are represented by the feature spaces. The linear scaling properties of the SMA fraction estimates from meter to kilometer scales also facilitate the vicarious validation of ET estimates using multiple resolutions of imagery

Spectral Characteristics of the Dynamic World Land Cover Classification

The Dynamic World product is a discrete land cover classification of Sentinel 2 reflectance imagery that is global in extent, retrospective to 2015, and updated continuously in near real time. The classifier is trained on a stratified random sample of 20,000 hand-labeled 5 × 5 km Sentinel 2 tiles spanning 14 biomes globally. Since the training data are based on visual interpretation of image composites by both expert and non-expert annotators, without explicit spectral properties specified in the class definitions, the spectral characteristics of the classes are not obvious. The objective of this study is to quantify the physical distinctions among the land cover classes by characterizing the spectral properties of the range of reflectance present within each of the Dynamic World classes over a variety of landscapes. This is achieved by comparing both the eight-class probability feature space (excluding snow) and the maximum probability class assignment (label) distributions to continuous land cover fraction estimates derived from a globally standardized spectral mixture model. Standardized substrate, vegetation, and dark (SVD) endmembers are used to unmix nine Sentinel 2 reflectance tiles from nine spectral diversity hotspots for comparison between the SVD land cover fraction continua and the Dynamic World class probability continua and class assignments. The variance partition for the class probability feature spaces indicates that eight of these nine hotspots are effectively five-dimensional to 95% of variance. Class probability feature spaces of the hotspots all show a tetrahedral form with probability continua spanning multiple classes. Comparison of SVD land cover fraction distributions with maximum probability class assignments (labels) and probability feature space distributions reveal a clear distinction between (1) physically and spectrally heterogeneous biomes characterized by continuous gradations in vegetation density, substrate albedo, and structural shadow fractions, and (2) more homogeneous biomes characterized by closed canopy vegetation (forest) or negligible vegetation cover (e.g., desert, water). Due to the ubiquity of spectrally heterogeneous biomes worldwide, the class probability feature space adds considerable value to the Dynamic World maximum probability class labels by offering users the opportunity to depict inherently gradational heterogeneous landscapes otherwise not generally offered with other discrete thematic classifications

Spatiotemporal Characterization of Mangrove Phenology and Disturbance Response: The Bangladesh Sundarban

This work presents a spatiotemporal analysis of the phenology and disturbance response in the Sundarban mangrove forest on the Ganges-Brahmaputra Delta in Bangladesh. The methodological approach is based on an Empirical Orthogonal Function (EOF) analysis of the new Harmonized Landsat Sentinel-2 (HLS) BRDF and atmospherically corrected reflectance time series, preceded by a Robust Principal Component Analysis (RPCA) separation of Low Rank and Sparse components of the image time series. Low Rank components are spatially and temporally pervasive while Sparse components are transient and localized. The RPCA clearly separates subtle spatial variations in the annual cycle of monsoon-modulated greening and senescence of the mangrove forest from the spatiotemporally complex agricultural phenology surrounding the Sundarban. A 3 endmember temporal mixture model maps spatially coherent differences in the 2018 greening-senescence cycle of the mangrove which are both concordant and discordant with existing species composition maps. The discordant patterns suggest a phenological response to environmental factors like surface hydrology. On decadal time scales, a standard EOF analysis of vegetation fraction maps from annual post-monsoon Landsat imagery is sufficient to isolate locations of shoreline advance and retreat related to changes in sedimentation and erosion, as well as cyclone-induced defoliation and recovery

Decades of urban growth and development on the Asian megadeltas

The current and ongoing expansion of urban areas worldwide represents the largest mass migration in human history. It is well known that the world's coastal zones are associated with large and growing concentrations of population, urban development and economic activity. Among coastal environments, deltas have long been recognized for both benefits and hazards. This is particularly true on the Asian megadeltas, where the majority of the world's deltaic populations reside. Current trends in urban migration, combined with demographic momentum suggest that the already large populations on the Asian megadeltas will continue to grow. In this study, we combine recently released gridded population density (circa 2010) with a newly developed night light change product (1992 to 2012) and a digital elevation model to quantify the spatial distribution of population and development on the nine Asian megadeltas. Bivariate distributions of population as functions of elevation and coastal proximity quantify potential exposure of deltaic populations to flood and coastal hazards. Comparison of these distributions for the Asian megadeltas show very different patterns of habitation with peak population elevations ranging from 2 to 11 m above sea level over a wide range of coastal proximities. Over all nine megadeltas, over 174 million people reside below a peak population elevation of 7 m. Changes in the spatial extent of anthropogenic night light from 1992 to 2012 show widely varying extents and changes of lighted urban development. All of the deltas except the Indus show the greatest increases in night light brightness occurring at elevations < 10 m. At global and continental scales, growth of settlements of all sizes takes the form of evolving spatial networks of development. Spatial networks of lighted urban development in Asia show power law scaling properties consistent with other continents, but much higher rates of growth. The three largest networks of development in China all occur on deltas and adjacent lowlands, and are growing faster than the rest of the urban network in China. Since 2000, the Huanghe Delta + North China Plain urban network has surpassed the Japanese urban network in size and may soon connect with the Changjiang Delta + Yangtze River urban network to form the largest conurbation in Asia